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1.   Introduction 

Populations  of  cells  that  make  up  organ  tissue  grow  and  contract  in  a  manner 
that  is  roughly  analogous  to  the  fluctuations  of  other  natural  populations.  Since 
organs  are  of  bounded  size  their  growth  is  not  entirely  uncontrolled  and 
exponential  overall,  or  otherwise  a  human  or  rat  liver  would  either  assume  a 
totally  outlandishly  large  size,  or  else  shrink  to  zero.  Neither  such  alternative  is 
seen  in  nature,  although  organ  sizes  do  vary  among  otherwise  comparable 
members  of  the  same  species.  And  organs  of  mature  hosts  can  change  in  size  as  a 
result  of  disease,  toxic  insult,  or  an  operation  such  as  partial  hepatectomy,  from 
which  a  liver  can  recover  again  to  normal  size  and  function. 


A  traditional  approach  to  modeling  population  (=  organ)  size  restriction  to  an 
observed  "normal"  level  is  to  postulate  a  carrying  capacity:  effectively  a  limit  on 
the  physical  region  that  can  be  occupied  by  the  population.  This  formulation 
apparently  goes  back  to  Verhulst  (1836);  see  Murray  (1989)  for  recent  discussion. 
In  the  organ  situation  this  might  correspond  to  a  space  of  approximately 
pre-ordained  dimension  that,  say,  liver  cells  in  liver  tissue  cannot  exceed  in  the 
body  of  a  mature  human  male.  The  space  can  be  taken  as  given,  introduced  into 
other  models  as  a  parameter,  and  in  particular  cases  estimated  from  data.  It 
would  be  the  maximum  size  of  the  liver  compartment  in  a  PB-PK  model,  for 
example. 

There  is  another  alternative  to  the  above  approach  that  depends  upon 
recognition  and  measurement  of  the  presence  of  various  biological  agents  called 
growth  factors.  There  are  a  number  of  such  factors  that  both  stimulate  (positively) 
and  inhibit  (negatively)  cell  population  growth.  Growth  factors  are  discussed  by 
Alberts  et  ah  (1994).  Aaronson  (1991)  provides  an  overview  of  growth  factors  in 
cancer;  see  also  Rubin,  Bottaro,  and  Aaronson  (1993).  The  purpose  of  this  note  is 
to  provide  a  very  simple  model  for  a  cell  population  that  grows  under  the  control 
of  positive  and  negative  growth  factors.  It  will  be  seen  that  such  a  model  can 
result  in  logistic  growth  without  the  necessity  of  postulating  a  physical  carrying 
capacity.  An  effective  carrying  capacity  appears  as  a  result  of  presumed  growth 
factor  interaction  with  cells. 

2.  Model  for  a  Cell  Population  Under  Growth  Factor  Control 

Suppose  a  population  of  cells  is  of  size  C{t)  at  time  t.  Its  individual  cell  growth 
or  birth  rate  is  Ao,  and  its  death  rate  is  jjq,  so  its  net  growth  rate,  Xq  -  [Uq,  governs 
the  manner  and  rate  of  growth.  Starting  with  C(0)  members,  and  left  alone,  the 
population  would  grow  roughly  like  C(t)  ~  C(0)e^~^°'  ,  which  means  either  to 


a  large  size  (Ao-^0  >  0),  or  to  zero  (Ao~  fJQ<0).  Clearly  such  unrestricted 
behavior  is  inappropriate  for  describing  a  population  of  cells  that  constitutes  an 
entire  organ,  although  essentially  such  a  model  has  been  used  to  describe  growth 
of  tumors  within  an  organ;  see  Tan  (1991)  for  an  overview;  in  particular  the  work 
of  Moolgavkar  and  co-authors,  cited  in  Tan  (1991),  is  relevant. 

Now  introduce  a  quantity  a(t)  of  a  positive  growth  factor  into  the  vicinity  of  the 
cell  population.  The  amount  present,  a(t),  changes  cell  birth  rate  to  Xq  +  X\ct(i), 
where  we  take  X\  >  0.  Also  introduce  a  quantity  /?(£)  of  negative  growth  factor;  it 
changes  cell  death  rate  to  jjq  +  Mpit),  ji\  >  0.  Then  changing  levels  of  a{t)  and  /X0 
can  certainly  alter  the  properties  of  the  cell  population,  from  net  growth  to  net 
decline,  depending  upon  values  of  exit)  and  p(t). 

Assume  that  the  productions  of  both  a(t)  and  j5(t)  are  regulated  by  cell 
activity  in  such  a  way  that 

^  =  p„C(t) -*„«(»)  (2.1) 

and 

That  is,  both  are  stimulated  to  increase  by  the  number  of  cells  present,  and  to 
diminish  in  proportion  to  their  own  concentration,  possibly  being  removed  from 
the  cell  site  (organ)  by  blood  flow  or  metabolism  or  other  biological  processes.  Of 
course  the  above  equations  are  prime  candidates  for  replacement  by  others  that 
more  accurately  depict  the  true  interactions. 

In  the  presence  of  a  and  j3-factors  the  cells  in  the  organ  grow  and  decline 
according  to 

£21  =  [A0  +  A,«(0]C(0  -  [no  +  MiflOFM.  <2.3) 


So  (2.1),  (2.2),  (2.3)  form  a  system  of  three  non-linear  differential  equations.  No 
explicit  solution  seems  immediately  available,  unless  we  make  the  quasi-static  or 
quasi-steady-state  assumption  (QSSA);  see  Strogatz  (1994)  for  its  invocation  so  as 
to  solve  a  non-linear  dynamics  problem  along  with  some  historical  references, 
and  Segel  and  Slemrod  (1989)  for  a  careful  discussion  of  this  approximation. 
Namely  assume  that  a(t)  and  p(t)  are  able  to  adapt  very  quickly  to  any  current 
value  of  C(t)  to  always  reach  a  "temporary  steady  state": 

fj  =  0  =  PaC(t)-Saa(t)  (2.4) 

^  =  0  =  ppC(t)-5p(3(t)-  (2-5) 

Adopt  the  approximation  as  true,  so  solve  (2.4)  and  (2.5)  for  ait)  and  p(t): 

oc(t)  =  (pa/Sa)C(t)  (2.6) 

and 

P(t)  =  (pp/8p)c(t).  (2.7) 

Let  us  call  (pa/  &a)  and  (pp/  8p)  the  prevalences  of  the  a  and  p  factors  respectively. 
Insert  these  into  (2.3)  and  for  convenience,  put  X{  =  A|(pa/<5a),  fi[  =  P-iippl^p), 


to  obtain 


^1  =  [^  +  A{C(0]C(0  -  [no  +  /iiC(r)]C(0 


=  (Ao-^)CW-(^-Ai)C2(0  (2.8) 

=  (Ao-/zo)C(0 


This  conforms  exactly  to  the  original  logistic  equation  if  the  ordinary  net  growth 
rate,  A  =  Xq  -  juq,  is  positive,  as  is  the  effective  carrying  capacity 


K  =  ^^£.  (2.9) 

Under  the  above  conditions  and  starting  from  C(0)  >  0,  the  population  attains  the 
long-run  steady-state  value 

CH  =  K  =  —( i?-^ -.  (2.10) 

lh\Pfil8p)-HPal8a) 

The  above  version  of  carrying  capacity  makes  intuitive  sense  in  that 

(a)  it  increases  with  net  population  growth  rate,  Ao  -  jjo; 

(b)  it  decreases  with  increased  prevalence  of  negative  growth  factor,  (pp/  8p), 
and  with  decreased  prevalence  of  positive  growth  factor,  ipa/  8q); 

(c)  the  inhibition  effect  of  negative  growth  factor,  \i\  =  Hiipp/dp),  must 
exceed  the  stimulative  effect  of  the  positive  growth  factor, 

M  =  hip al &ob- 

If  any  of  the  above  conditions  are  violated  the  population  development 
becomes  radically  different,  but  can  also  be  interesting. 

The  time-dependent  population  size  is  seen  to  be  of  the  familiar  logistic 
growth  form 


KC(0)eAt 

K-c(o)+c(oy 


CM--     ^hzTM  (2-lD 


with  K  as  in  (2.10),  A  =  Ao  -  /io  >  0,  and  0  <  C(0)  <  K. 

Note  that  the  formula  has  biological  meaning  even  if  C(0)  >  K,  and  also  if 
K  <  0:  suppose  that  A  =  Ao  -  fJQ  >  0  but  A{  >  ji{;  then  write  K'  =  -K  >  0  to  get 

^l  =  AC(t)[l  +  C(t)/K']f  (2.12) 

the  solution  to  which  is 


K'C(0)eM 

K'+c{o)-c(oy 


**)=„  Jzz«**  <2-13> 


if  t<-\n 
A 


1  + 


v 


C(0) 


1 
;  it  explodes  when  t  =  —  In 


1  + 


C(0)y 


;  this  might  plausibly 


model  an  especially  malignant  tumor  growth.  Finally  when  A  =  Ao  -  fjQ  <  0  and 
K'  =  -K  >  0  we  simply  get  (2.11)  with  -A  replacing  A,  once  again  a  logistic  model, 
but  now  one  that  decreases  as  t  increases. 

3.   Stochastic  Models 

The  above  model  can  be  "made  stochastic"  in  various  ways.  One  is  to  re-state 
the  growth  factor  and  cell-growth  equations  as  a  system  of  three  non-linear 
Ito-type  stochastic  differential  equations.  Analytical  solutions  are  not  likely  to  be 
available,  but  some  asymptotics  might  well  produce  explicit  results. 

Another  approach  is  to  assume  that  a  stochastic  version  of  C(t),  namely  C(t),  is 
a  birth-and-death  process  with  transition  rates  copied  from  the  right-hand  side  of 
(2.8).  For  example,  C(t)  evolves  over  state  space  (0, 1, 2, ...)  according  to 

P{C(t  +  dt )  =  C{t)  +  l\C(t)  =  i)  =  (Xq  -  ii0)idt  +  o(dt)  (3.1,a) 

P{C{t  +  dt)  =  C{t)  -  l\C(t)  =  i)  =  (/zj  -  X^fdt  +  o(dt);  (3.1,b) 

these  holding  when  A  =  Ao  -  ^o  >  0/  Hi  -  -M  >  0;  otherwise  modification  is  needed. 
Although  E\C(t )|C(0)  =  C(0)1  *  C(t )  because  of  the  non-linearity  in  the  generator, 

(3.1),  it  can  be  of  interest  to  study  the  above  stochastic  version's  transient 
properties,  such  as  first-passage  times  from  low  states  (population  size)  to  high, 
or  the  reverse,  e.g.  to  C(f#)  =  0  when  the  population  dies  out. 
The  described  approach  essentially  minimizes  attention  to  the  stochastics  of  the 
growth  factors  and  ignores  non-linearity,  hence  is  a  prime  candidate  for  an 


upgraded  treatment.  Nevertheless  it  is  appealing  for  its  simplicity  and  easy 
availability,  and  is  offered  as  an  interim  approach. 

4.   Summary 

It  is  shown  that  classical  logistic  growth  can  be  induced  in  a  non-traditional 
manner  by  hypothesized  action  of  growth  factors,  rather  than  by  action  of  a 
physical  carrying  capacity  (although  the  latter  may  operate  as  well).  Modification 
of  the  effective  carrying  capacity  to  be  negative  has  biological  interpretability. 
The  resulting  models  may  perhaps  find  a  use  in  cell  proliferation  and  cancer 
modeling. 
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Dept.  of  Operations  Research 
Stanford  University 
Stanford,  CA    94305 

32.  Prof.  Linda  V.  Green 

Graduate  School  of  Business 
Columbia  University 

New  York,  NY    10027 


33.      Prof.  J.  Michael  Harrison 

Graduate  School  of  Business 
Stanford  University 
Stanford,  CA   94305-5015 


34.  Dr.  D.  C.  Hoaglin 

Department  of  Statistics 

Harvard  University 
1  Oxford  Street 
Cambridge,  MA   02138 

35.  Dr.  David  G.  Hoel 

Professor  of  Biometry  and  Epidemiology 
Medical  University  of  South  Carolina 
171  Ashley  Ave. 

Charleston,  SC   29425-0002 


36.  Prof.  D.  L.  Iglehart 

Dept.  of  Operations  Research 
Stanford  University 
Stanford,  CA    94305-5015 

37.  Institute  for  Defense  Analysis 
1800  North  Beauregard 
Alexandria,  VA   22311 


38.      Dr.  Robert  C.  Jackson 

Vice  President,  Research  and  Development 
Agouron  Pharmaceuticals,  Inc. 
3565  General  Atomics  Court 
San  Diego,  CA   92121-1121 


39.      Prof.  J.  B.  Kadane 

Dept.  of  Statistics 
Carnegie-Mellon  University 
Pittsburgh,  PA    15213 


40.      Dr.  F.  P.  Kelly 

Statistics  Laboratory 
16  Mill  Lane 
Cambridge,  ENGLAND 
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41.      Dr.  Jon  Kettenring 

Bellcore 

445  South  Street 

Morris  Township,  NJ   07962-1910 


42.      Dr.  Ralph  Kodell 

Chief,  Biometry  Branch 
Biometry  and  Risk  Assessment 
National  Center  for  Toxicology  Research 
3900  NCTR  Drive 
Jefferson,  AR   72079 


43.  Mr.  Koh  Peng  Kong 

OA  Branch,  DSO 
Ministry  of  Defense 

Blk  29  Middlesex  Road 
SINGAPORE  1024 

44.  Prof.  Guy  Latouche 

University  Libre  Bruxelles 
C.P.  212,  Blvd.  De  Triomphe 
Bruxelles    B-1050,  BELGIUM 


45.  Dr.  A.  J.  Lawrance ... 

Dept.  of  Mathematics 
University  of  Birmingham 

P.O.  Box  363 

Birmingham    Bl 5  2TT,  ENGLAND 

46.  Prof.  M.  Leadbetter 

Department  of  Statistics 
University  of  North  Carolina 
Chapel  Hill,  NC   27514 


47.      Prof.  J.  Lehoczky 

Department  of  Statistics 
Carnegie-Mellon  University 
Pittsburgh,  PA    15213 


48.  Dr.  Georg  Luebeck 

Fred  Hutchinson  Cancer  Research  Center 
1124  Columbia 

MP-665 

Seattle,  WA   98014 

49.  Dr.  Colin  Mallows 

AT&T  Bell  Telephone  Laboratories 

600  Mountain  Avenue 
Murray  Hill,  NJ   07974 
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50.  Prof.  R.  Douglas  Martin 

Department  of  Statistics,  GN-22 
University  of  Washington 
Seattle,  WA   98195 

51.  Dr.  Sari  Mazumdar 

Biostati  sties  Dept. 

University  of  Pittsburgh 
Graduate  School  of  Public  Health 
Pittsburgh,  PA    15261 

52.  Dr.  James  McKenna 

Bell  Communications  Research 
445  South  Street 

Morristown,  NJ   07960-1910 

53.  Dr.  Ramit  Mehr-Grossman 

Theoretical  Biology  and  Biophysics 
Theoretical  Division 

Mail  Stop  K710 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM   87545 

54.  Prof.  Carl  N.  Morris 

Statistics  Department 

Harvard  University 
1  Oxford  Street 
Cambridge,  MA    02138 

55.  Dr.  John  A.  Morrison 

AT&T  Bell  Telephone  Laboratories 
600  Mountain  Avenue 

Murray  Hill,  NJ    07974 

56.  Prof.  F.  W.  Mosteller 

Department  of  Statistics 
Harvard  University 

1  Oxford  Street 
Cambridge,  MA    02138 


57.  Dr.  John  Orav 

Biostatistics  Department 
Harvard  School  of  Public  Health 
677  Huntington  Avenue 

Boston,  MA   02115 

58.  Dr.  Alan  Perelson 

Theoretical  Biology  and  Biophysics 
Theoretical  Division 

Mail  Stop  K710 

Los  Alamos  National  Laboratory 

Los  Alamos,  NM    87545 
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59.      Dr.  Jim  Petty 

National  Biological  Survey 
4200  New  Haven  Road 
Columbia,  MO   65201 


60.      Dr.  Lorenz  Rhomberg 

Harvard  Center  for  Risk  Analysis 
Harvard  University 
Cambridge,  MA    02138 


61.  Dr.  Rhonda  Righter 

Dept.  of  Decision  &  Info.  Sciences 

Santa  Clara  University 
Santa  Clara,  CA   95118 

62.  Dr.  John  E.  Rolph 

Information  and  Operations  Management 
Univ.  of  Southern  California 

School  of  Business  Administration 
Los  Angeles,  CA   90089-1421 


63.      Prof.  M.  Rosenblatt 

Department  of  Mathematics 
University  of  California,  San  Diego 
LaJolla,CA   92093 


64.      Prof.  Frank  Samaniego  .. 
Statistics  Department 
University  of  California 
Davis,  CA   95616 


65.  Prof.  G.  A.  F.  Seber 

Dept.  of  Statistics 

Univ.  of  Auckland 
Private  Bag  92019 
Auckland,  NEW  ZEALAND 

66.  Prof.  G.  Shantikumar 

The  Management  Science  Group 
School  of  Business  Administration 
University  of  California 
Berkeley,  CA   94720 

67.  Prof.  N.  D.  Singpurwalla 

George  Washington  University 
Washington,  DC   20052 

68.  Prof.  H.  Solomon 

Department  of  Statistics 
Sequoia  Hall 

Stanford  University 
Stanford,  CA    94305 
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69.  Dr.  Andrew  Solow 

Woods  Hole  Oceanographic  Institute 
Woods  Hole,  MA   02543 

70.  Prof.  W.  Stuetzle 

Department  of  Statistics 
University  of  Washington 

Seattle,  WA    98195 


71.  Prof.  J.  R.  Thompson 

Dept.  of  Mathematical  Science 
Rice  University 

Houston,  TX    77001 

72.  Prof.  Steven  K.  Thompson 

Statistics  Dept. 
Pennsylvania  State  Univ. 

326  Classroom  Bldg. 
University  Park,  PA    16802-2111 


73.      Prof.  J.  W.  Tukey 

Statistics  Dept.,  Fine  Hall 
Princeton  University 
Princeton,  NJ   08540 


74.      Dr.  D.  Vere-Jones 

Dept.  of  Math 

Victoria  Univ.  of  Wellington 

P.O.  Box  196 

Wellington,  NEW  ZEALAND 


75.      Prof.  David  L.  Wallace . 
Statistics  Dept. 
University  of  Chicago 
5734  S.  University  Ave. 
Chicago,  IL  60637 


76.  Dr.  Ed  Wegman 

George  Mason  University 

Fairfax,  VA    22030 

77.  Dr.  L.  Wein 

Operations  Research  Center,  Rm  E40-164 
Massachusetts  Institute  of  Technology 
Cambridge,  MA    02139 

78.  Dr.  Alan  Weiss 

Rm2C-118 

AT&T  Bell  Laboratories 
600  Mountain  Avenue 
Murray  Hill,  NJ    07974-2040 
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79.  Prof.  Roy  Welsch 

Sloan  School 

M.I.T. 

Cambridge,  MA    02139 

80.  Dr.  Raymond  S.H.  Yang 

Colorado  State  University 

College  of  Veterinary  Medicine  and  Biomedical  Sciences 
Dept.  of  Environmental  Health 
Fort  Collins,  CO   80523 
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